Wavelets and Quantum Algebras 



A. Ludu^, M. Greiner* and J. P. Draayer^ 

^Department of Physics and Astronomy, Louisiana State University, 
Baton Rouge, LA 70803-4001, U.S.A. 
^Institut fur Theoretische Physik, Technische Universitat, 
D-01062 Dresden, Germany 

Abstract: 



Wavelets, known to be useful in non-linear multi-scale processes and in multi-resolution 
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Haar and B-wavelets are worked out in detail. 
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1 Introduction 



Interest in wavelet theory and its applications in multi-resolution analysis [1] has grown over 
the last decade, involving new and linking disparate fields of research from pure mathe- 
matics and physics to down-to-earth signal engineering. It is being widely applied in signal 
processing and data compression [2-4] , pattern recognition [5] , statistical physics and turbu- 
lence [6], jet dynamics [7], field theory [8,9], solid state physics [10-12], quantum mechanical 
applications [13-16], non-linear dynamics [6,7,17], soliton theory [18,19], etc. 

The rise in interest in wavelet theory is motivated by the fact that Fourier analysis is 
ineffective when dealing with non-linear models and localize sharp features [6-18]. A central 
property of wavelets is their ability to expand and analyze functions with respect to a set of 
self-similar localized basis functions (scaling functions and wavelets) and to processes them 
locally without affecting the scale. The wavelet method is recursive and therefore ideal for 
computational applications. Moreover, the scaling functions and the corresponding wavelets 
are very well localized both in the time and frequency domains, and hence wavelets are ideal 
for an analysis of phenomena where different space/time scales occur. They also provide 
mathematical representations that can handle both analytical and numerical difficulties due 
to singular phenomena [6,7]. Like Fourier analysis, wavelet theory uses basis functions with 
different characteristic scales. However, Fourier analysis has the advantage of being build 
upon a simple and solid mathematical foundation [20]. The challenge is to create such a 
platform for wavelet theory. 

This paper is a first step in that direction. It focuses on obtaining a closed algebraic 
method for the construction of the scaling and wavelet functions. Results for the Haar and 
B-wavelets are given. The results show that in addition to providing a multi-resolution 
basis of self-similar functions, wavelets display a definite non-linear symmetry which can be 
associated with a q-deformation of the Fourier series generating algebra. In a second paper 
we will show that wavelets also have variational properties: all multi-scale equations follow 
from Hamilton's equation of an infinite-dimensional Hamiltonian system. 

The two-scale equation, which is the central object of wavelet theory, is also shown 
to have an algebraic structure that is associated with a quantum algebra [21]. Applications 



2 



of this result can be realized using supercomputers and efficient numerical schemes which, in 
turn, can be used to model complex continuous systems or genuine discrete systems defined 
on lattices. The symmetries related to complete solvability [22] or to discrete-continuous 
transitions [23], and their basic tools, the finite difference equations, are new and key elements 
in any analysis of discrete systems. On uniform lattices the symmetry algebras of partial 
differential equations is left unchanged by the discretization [24] . In the case of non-uniform 
lattices one needs to introduce generalized symmetries, like quantum algebras [25]. This is 
an example of when the q-deformation of some initial symmetries can play an important 
role. An outcome of the present paper is an extension of traditional Fourier analysis towards 
wavelet expansions by means of q-deformation. 

The finite-difference and scaling operators involved in the dilation equation are alge- 
braically closed with respect to certain non-linear commutation relations. This symmetry 
is seen most simply when the operators are realized in terms of q-deformed derivatives [26]. 
By expressing the translation and dilation operators as q-deformed derivatives, the scaling 
and difference operators can be mapped onto the generators of a non-linear algebra and the 
action of the q-deformed derivatives is extended to non-different iable functions. Moreover, 
such non-linear algebras can be un-deformed to the Fourier series generating algebra or can 
be mapped onto the su q {2) algebra. 

After introducing notation and basic definitions in Section 2, we present in Section 3 
the deformation of a Fourier algebra into a scaling function generating algebra, and identify 
a spectral problem with the Haar dilation equation. In Section 4, we introduce the defini- 
tion of the most general scaling function and wavelet algebra and prove that the algebraic 
formulation of the dilation equation is unique. A duality relation is obtained: for any scaling 
function there is a non-linear algebra; and for special algebras one can find a corresponding 
wavelet structure. In Section 5 we give examples of scaling function algebras for the Haar and 
the B-wavetes. We use new limiting procedures to solve the q-difference and finite-difference 
equations. In Section 6 we use a variant and extension of the deforming functional technique 
[28,29,31] to obtain a mapping of the generators of non-linear algebras to those of su q (2). 
Concluding comments, further extensions, and other remarks are given in Section 7. 
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2 Basic elements 



An algebraic structure can be built for any scaling function system and wavelet basis. In 
order to realize this construction, three basic building blocks are needed: the algebraic 
structure of the Fourier system (starting point), the wavelet theory (final objective), and 
q-deformation (intermediate tool). 

2.1 Fourier algebras 

In order to demonstrate our approach, consider algebras which have Fourier series as the basis 
of their representation spaces [20]. In the following we shall denote d k f/dx k = d^f = f k . 
The trigonometric (Fourier) system {\k >= e lkx }kez diagonalizes all translation invariant 
operators acting on L 2 ([0, 2ir]). We introduce three generators within a differential realization 
Jo = —id, J± = e ±%x {— id) p , satisfying the commutation relations 

[J ,J±] = ±J±, [J + ,J-\ = l-{l-2idf. (1) 

For p — 0, eqs.(l) describe an algebra, denoted Tq, that is isomorphic to an analytical 
prolongation, e(2,C), of the Euclidean algebra e(2,R) of rigid motions in the plane. The 
algebra e(2,R), generated by P x = P^, P y = Pj and R = —R\ with the commutators 
[R, P x , y ] = ±P y ,x, [Px, Py] = 0, is realized through the mapping 

Jo = iR, J± — Px i iPy, 

onto the Tq algebra. The unitary irreducible representations (unirreps) of are based on the 
eigenvectors of the self-adjoint operator J , Jo\n >= n\n >. For any two distinct eigenstates, 
| ra > and \n >, by using the first commutator in eq.(l), we obtain n — n + 1 and hence 
the spectrum of Jo is unbounded, discrete and consists in equidistant eigenvalues. Hence 
the space of representations is generated by the Fourier system, and the other generators 
J ± |n >= \n ± 1 > act like ladder operators on the \n > states, increasing/decreasing the 
scale by unity. For p — 1 eqs.(l) describe another Lie algebra, T\, isomorphic with the 
symplectic algebra sp(2,R) ~ su(l, 1) ~ so(2, 1). This differential realization has the same 
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representation space as JF but it is not irreducible. The action of the generators of T\ is 
similar to that of the generators of JF for n ^ 0. Unlike the JF case, however, we have that 
J± \ =Fl >= |0 > and the sub-spaces {e mx } n€N and {e~ mx } n€N are invariant subspaces of T\. 
The JF 01 examples suggest the possibility of other constructions in terms of the operator d 
and the complex exponential functions of different scales. 



2.2 Scaling functions and wavelets 



Wavelets are able to reconstruct a signal through regular sampling, namely, by analyzing the 
signal at different scales (which increase/decrease exponentially) with the step size between 
each scale being the same. Different from the ordinary Fourier transform, which reproduces a 
function as a superposition of complex exponentials, or from the windowed Fourier transform, 
which introduces a scale into the analysis of signals, multi-resolution analysis processes the 
signal locally, using the appropriate local scale. Basically, wavelets are constructed with a 
pair of operators: the dilation (scaling) and finite-difference (combinations of translations) 
operators, acting on L 2 (R) and defined, respectively, by T a f(x) = f(x + a), D^fi^x) = 
2 l3 f(2 /3 x), with a, f3 arbitrary real numbers. They are invertible, unitary, and fulfill T a D@ = 
D P T 2Pa^ T a T f3 = T a +I 3 ? D a D p = D a+f3 _ The formal Tay i or series representation of these 

operators are T a = e ad and D 13 = 2^e^ ln2xd . On the space of compact supported or rapidly 
decreasing functions, any holomorphic function f(T) is locally polynomial. Indeed, if / 
is holomorphic, and its action is taken on the compact supported function subspace of 
$ G L 2 (R), then the action of f(T) = J2k&zCkT k reduces to that of a Laurent polynomial 
by keeping only a finite number of terms in the sum. 

The wavelet system is given by a set of scaled and translated copies of a pair of 
functions: the scaling function $ and the mother wavelet The basic fact about wavelets 
is that both these fundamental functions are finite linear combinations of reflecting the 
self-similar character of the wavelet system. The defining equation for the scaling function 
(the dilation equation) is a linear finite-difference equation, including a scale change (q- 
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difference [27]) 

n 

$(x) =DJ2 C k T- k <Z>(x) = Dg(T)$(x), (2) 

k=0 

where the RHS sum is a polynomial in T, g(T). Eq.(2) is a fixed-point equation and conse- 
quently it has only one unique solution [1-7,11,13]. The scaling function $(x), as a solution 
of eq.(2), is required to have two properties [1-5]: 

1. f R $(x)dx = 1, the average value property; 

2. < $„, $ >= j R &(x + n)Q(x)dx = S n fi for any n E Z, the orthogonality condition. 
These conditions introduce restrictions on the coefficients C k in eq.(2) [1-4], 

n n 

^ Ck = 2, ^ C k C k+ 2i = S i, I G Z. (3) 

k=0 k=0 

The conditions expressed through eq.(3) yield a pattern of L 2 (R) as a chain of subspaces 
Vj C Vj+i, each one being generated by all the translations of D j &, j e Z. By repeated 
application of -D and T on <£>, D^T n ^{x) = $(2 J x + n) = one obtains a non-orthogonal 
basis in L 2 (R). The action of the operator —DT~ x g(—T~ 1 ) on $, A being a unique odd 
integer, provides the mother wavelet function: 

V(x) = -T^ZM-T- 1 )^) = -T- A ^(-l) fc C_ fe+1 $(2 3 ; - fc). (4) 

The wavelet ty(x) has the property that {^j, n }j,nez is an orthonormal basis of L 2 (R) [1-3] 
and L 2 (R) is a direct sum of the orthogonal subspaces Wj (the orthogonal complements 
of Vj), each of them generated by all possible translations of = with integral 

m: L 2 (R) = ®j e zWj, [1-5]. Hence translated and dilated copies of the mother wavelet, 
^{QPx + k), generate a true orthonormal basis. Recursively applications of the eqs.(2,4) 
relate all the scaling functions and wavelets. 

The simplest example is provided by the Haar wavelet, defined by the scaling function 
^Haar(x) — 1 if \x — 1/2 1 < 1/2 and otherwise. In this case we have D(l + T~ l )§ H aar = 
&Haar and gHaar(T) = 1 + T -1 (C = C-\ = 1 and the rest in eq.(2)). The corresponding 
Haar wavelet is ^ Haar = D(l — T~ 1 )& Haar , i.e. A = 1 in eq.(4). 

There are similarities and differences between the Fourier and the wavelet approaches. 
From an algebraic point of view, in both approaches, there are eigenfunction equations, which 
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keep the scale constant, and the ladder operators, which change the scale. The dilation 
equation, eq.(2), has its analog in the Fourier formalism, though it is a fixed-point equation 
and has only one solution. The mother wavelet, eq.(4), has no analog in Fourier analysis. 
Each Fourier eigenfunction carries one scale; a wavelet function ^/(x) involves two scales, e.g. 
$o,fc an d $2,fc- Because of their localization, wavelets possess a degree of freedom beyond 
that of a Fourier system; namely, the width of the support of $. In the present approach 
this degree of freedom is associated with the deformation parameter, q. 



2.3 Quantum deformations 



A possibility for constructing scaling functions/wavelets within an algebraic approach is to 
deform the Fourier algebraic structure into a non-linear system. In general, the scaling func- 
tions and wavelets are not differentiable. This suggests the use of finite-difference operators 
instead of derivatives. Finite-difference operators are closed with respect to commutation but 
only within non-linear algebraic constructions. Consequently, it is natural to use q-deformed 
derivatives to find a foundation for wavelet theory in g-deformed algebras. 

Quantum algebras refer to some specific deformations of Lie algebras, to which they 
reduce when the deformation parameter q is set equal to unity (for a recent monograph see 
[21] and references therein). The simplest example of a g-algebra is su g (2) whose Jordan- 
Schwinger realization is given in terms of g-bosonic operators [26]. The q-deformed algebras 
have been applied in various branches of physics like spin-chain models, non-commutative 
spaces, rotational spectra of deformed nuclei, Hamiltonian quantization, and dynamical sym- 
metry breaking. The basic element is the q-deformation of a certain object x, which can be 
a number, an operator, or a function: 

Ms = - — ^ x > ( 5 ) 

q-q 1 

where q = e s . The Taylor expansion of T and D are related to the q-deformation of the deriva- 
tive operator, according to the definition of the coordinate description of the q-deformed 
oscillator introduced in [26], 



Analogously, we introduce the operator 

r*0i, f (x) - f{2Sx) ~ f{2 ~ Sx) - DS ~ D ~ S f(x) (7) 

When s — > 0, [d] s — > d and [x£?] s in2 — > Eqs.(6,7) can be inverted and the occurrence of 
the q-deformed derivative and the translation/dilation operators becomes immediate 

T ±s = -(±% + ^- x +2] (8) 



2 V rj(s) r] 2 (s/2) 

D ±s 1 / , [Xd] a M 2 s/2 N 

"2l ± )] (sln2)% 2 ( S ln2/2) + ^ W 
where s £ N (or more general G -R) and 77(5) = eS _ 1 e _ s . The q-deformed algebra su q (2), 
generated by {Jo, J±}, is defined by the q-deformed version of the commutation relation 

[J + ,J_] = [2J ] S 

while the other two commutator relations remain undeformed, that is, they have the same 
form as for su{2). In addition to this traditional version of su q (2), several generalized forms 
have been introduced through two different prescriptions: 1) by deforming the commuta- 
tor [Jo, J±] by using some arbitrary function G(J ,q), [J , J+] = G(J ,q)J + and [J , J_] = 
— J_G(J ,g), independently proposed in [28], [29] and [30]; and 2) by deformations involv- 
ing all three commutation relations, using two functions G(J ,q) and F(J ,q) = [J + , J_], 
introduced in [31]. Unlike the former, for which the spectrum of Jo is linear, the latter is 
characterized by an exponential spectrum for Jo. Since in wavelet theory the domains of 
analysis are divided exponentially, rather than linearly with bands of equal widths, such 
algebras are used in the following analyses. 



3 Haar scaling function algebra, A s , 



The aim in this section is to obtain operators depending on D and T that form a non- 
linear algebra. This structure must provide an algebraic form for eqs.(2,4) and map onto the 
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Fourier algebra, JF i. Consider an operator depending on T and two real parameters s, a: 

x T Sa - T- Sa COS S7T . . 

Wo (s, a) = ^p-r , 10 

24(a) smhs 

where £(a) = sin ^ sin(^ L ) — 2icos(^ E ). The (s, a) parameters allow to approach d or 
combinations of T 

W (0,a) = — -9, (11) 
?(«) 

W (-,a) = — ^— — -Tt, (12) 

U V2' / 2sinh(l/2)f(a) ' v ; 

W (l,a)= , (13) 
2 smh(l)4(o;) 

w ^ a) - 2 SmHma y (14) 

In the limit s — > 0, W reduces to the normal derivative with respect to x, eq.(ll). In 
the limit s — |, we obtain a power of the translation operator T, eq.(12). Eq.(13) de- 
fines a q-deformation of unity, namely, Wo(l,0) = 4sin t h( - 1 - ) . Wo(s,a) can be a Hermitian or 
anti-Hermitian operator, W (l,a) = ±W (l,a)\ depending on a. The last limit, eq.(14), 
represents a finite-difference operator and is proportional to the q-derivative with respect to 
the a-deformation [2d] a , similar to eqs.(6-9). For a = 2k in eq.(12), a = k in eq.(13) and 
a = k/2 in eq.(14), the corresponding W (s,a) are Laurent polynomials in T, providing a 
direct connection to the wavelet operators. For some values of a, eqs. (12-14) give invertible 
operators with respect to T, T = T(Wq). 

In order to introduce the dilation operator D, we define 

W±(s) = -e^^D^e^^il + T ±s ). (15) 

In order to fulfill the hermiticity condition (W^) 1 " = W- on L 2 (R) these generators can be 
redefined in the form 

W±^W± = 2 ^2e T%x ~D Ts e T%x ~ I 1 + e % 4 y±s(2z ) 

By using the commutators relations between D and T and by performing an integration by 
parts, one has a direct check of the relation = W- 

< fi, W+h >= I fi(x)W + f 2 (x)dx =< W-h, h > . 
Jr 



From eqs.(10,15) one obtains the commutators 

[W (s, a), W + (s)\ = G(s, T)W + (s), (16) 

\W (s,a),W-(s)] = -W-(s)G(s,T), (17) 

\W + (s),W-(s)]=F(s,T), (18) 

with 

G(s, T) = W (s, a) ^ 

and 

F(s, T) = - ( (1 + e i^ s )^T 2Ss )(l + T~ s ) - (1 + gi-a+a-J^r-a--) (! + T *<\ _ ( 2 0) 
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Eqs. (10, 15, 16-20) describe a non-linear associative algebra, generated by T ± ,D ± , denoted 
A S) a- The eigenvalue problem for W± provides the algebraic form for the scaling equation. 
When Wq(s, a) is invertible with respect to T, the algebra has Wo, W± as generators which 
depend on two parameters (s, a) and is homomorphic with a special q-deformed algebra, 
namely, the two-color quasitriangular Hopf algebra A ± [31]. The spectrum of W is not 
equidistant because F{s, T) ^ const, x W . The Casimir operator of the algebra A s , a is 
given by [28-31] 

C = W_W + + H(W ) = W + W_ + if (W ) - F(W ), (21) 

where if is a real function, holomorphic in a neighborhood of and which must satisfy the 
functional equation 

H(t)-H(t-G(t)) = F(t). (22) 

with £ generic. The irreps of A s , a & re labeled by the eigenvalues W \a >= a\a >. The 
commutator in eq.(16) can be written in the form 

W+W Q = (W - G)W+. (23) 

If | a >t^ | a' > are eigenvectors of Wo, with corresponding eigenvalues a, a', we have from 
eq.(23) a non-linear recursion relation for all the eigenvalues 

a=a-G(a). (24) 
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The algebra A s , a can be mapped into su(l,l), the Lie algebra of a Fourier series, in the 
limit s — > 0, a — > 2: W± — > e ±i:r = J±, W — > — i<9 = J , where we take these limits in 
the order »4 S0iQ0 = lim a ^ ao lim s „» so A s , a - The continuous mapping A s , a ~ * ^0,2 — su(l, 1) is 
an algebraic morphism and one can consider wavelets as q-deformed generalizations of the 
Fourier series, in the above sense. 

As a first example, the Haar scaling function can be associated with a particular case 
of A s , a , namely, A\ t \. In the algebra A\ t i the eigenproblem for W_ provides the dilation 
equation for the Haar scaling function $(x) defined by eq.(2) with g(T) = 1 +T -1 . We have 
in this case G — W — 2W$ + 1, and W± = \D^ l {\ + T ±r ) and the commutators 

[w 0: w + ] = (Wo - ml + i)w+, 

[Wo, W-\ = -W-(W - Wl + 1), (25) 

[w+,wq = 1 -(t\Wo)-t 1 i\w )-t- 1 i 2 (Wo) + t- 1 (Wo)^, 

where T(Wq) is the solution of the equation 2Wq = T + T -1 , a pseudo-differential operator. 
The Casimir operator of ^4i i is a constant. Indeed, eq.(22) for the general Casimir operator 
is given by 

H{Wo)-H{2Wl-l) =F(W ), (26) 
and has, in terms of the T = T(W ) operator, the form 

X (T) - X (T 2 ) = \{t 2 - T 1 ' 2 - T- 1 ' 2 + T- 1 ) , (27) 

where xCO — H( T+ ^ ). The unique analytical solution for eq.(27) reads 

X (T) = i (^const. - T - T l/2 - T~ 1/2 ) , (28) 

which gives for the Casimir operator a constant. 

The spectrum of Wq consists of periodic functions satisfying the equation <3>(x + 1) + 
<fr(x — 1) = 2$ (x). This takes one back to Fourier analysis. The way to the Haar scaling 
function is to use the spectrum of W± since in A\ t \ the eigenproblem for W- gives the 
dilation equation for the Haar scaling function Once <3> is obtained from = $ 

the generator W carries <£> into an infinite sequence of functions Wq&. These functions 
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are mutually orthogonal and generate the space V . The action of W± on $ yields Wq"&. 
Consequently V is an invariant space for all the generators of the algebra A\^. 

The spectrum of W depends on the values of the parameters s, a. In the case of 
Ai,\ this relation is a = 2a 2 — 1. The corresponding representations are infinite-dimensional 
and the eigenvalues of W | 1,1 are a n = cosh(2 n a), n G Z for any a. The corresponding 
eigenfunctions have the form \a n >±= e ±2 " ax which results in an exponential spectrum for 
a n similar to the sequence of scales in wavelet theory. This is a self-similar spectrum with 
respect to a, like the sequence of scales in wavelet theory. A part of this spectrum is shown 
in Fig. 1. The action of W± is: 

1 + e ±2 " d 

W±\a n >±= |a n± i > ± . (29) 

Another basis for a representation of Ai : i is given by \k > = e 1 ^, k £ Z. We have the action 

W \k > = cos j jk > , W±\k > = l - (\ + e^ ^A; > , 

and on this basis (W r +) t = W_. There are two invariant spaces, \2 k > and \2~ k >, k G N, 
with W±\0 > = 1. 

The same procedure can be followed for any set of the parameters s, a. For example, 
if we choose the algebra ^.2,1/2 with W defined in eq.(12), we obtain the spectrum of W 
described by the recursion relation: 

a =a(a + Va 2 + 1 + \= \ (30) 

V a + V a 2 + 1 / 

which also gives a non-linear, unbounded, discrete representation of ^.2,1/2- 

The last thing to prove is the uniqueness of the two-scale equation in the A± t i scaling 
function generating algebra. The two-scale equation, eq.(2), in its algebraic form 2W^_$ = $ 
is not unique in Ai t i if there exists an operator, similar with that occuring in the dilation 
equation, which commutes with W_. This operator should contain higher powers in T and 
consequently is an element of U(Ai t i). 

Proposition 1: In U(Ai t i) there exists a unique operator X = Dx{T~ x ) such that 
\W-,X\ = 0. The function x(£) is integer, not polynomial, and unbounded, for generic 

Proof: We take for X a Laurent series x(T _1 ) = J2kez C k T~ k then the condition 
= results in a recursion relation for the coefficients C k : 

C2fc+1 = Ck — C2fc-1, C*2fe = Ck — C2k-2, (31) 
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for any k £ Z. Eqs.(31) have one trivial solution which reproduces W-, x(T 1 ) = C_i(l + 
T _1 ) and only one other solution, with C k = ±C_i or 0, uniquely defined (C 2 = C± 3 = 



-4 — ^5 — ^-6 — <-/7 — 0_g — U_n... — U, U_2 — <^-5 — — ^-9 — ^-10 — WO — ^11 — 



... = C_i and Ci = C 4 = C_ 7 = C 8 = C g = C_ 12 = C i2 = ... — C_i, etc.). The sequence of 
non-zero coefficients is infinite, and {C^} is not a Cauchy sequence, (q.e.d.) We note that 
the Ai,i algebra is also a Hopf algebra, [21,30,31], defined by eqs.(25) and by the coproduct, 
counit and antipode in the form: 



All its generators are primitive elements. 



4 General scaling and wavelet algebra 

This algebraic approach for the Haar scaling function can be generalized to yield a non-linear 
algebra for any scaling function, and conversely, to find the dilation equation and scaling 
function for certain algebras. The procedure starts with an algebra generated by the dilation 
and translation operators and constructs, within this algebra, the two-scale equation. The 
generators of the A s , a algebra can be generalized as 



with jo(T) and j(T) being arbitrary functions of T and s, holomorphic in a neighborghood 
of T = 1, with their dependence on s being such that in the limit s — > 0, jo(T) — > —id, 
j(T) — > 1. In this case the commutation relations, eqs. (16-20), become 



AT ± = T ± ®T ± , AD ± = D ± ®D ± , e^) = e(D ± ) = 1, 



W ^j (T), W ± ^j ± {D,T) = e^^D^e^^j^ 1 ), 



(32) 



[jo,j+] = Gj + , \j ,j-\ 



■j-G, [J + ,J-]=F 



(33) 



with G, F depending on T through jo,j, respectively, 




(34) 
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F = F(s, T) = j(s, e^^T* s )j( S , T _1 ) - j(s, T" 2 "' e^^)j(s, T). (35) 



Eqs. (33-35) define a non-linear algebra denoted Aj j as a generalization of A s , a - If the 
function j(T) is invertible with respect to T, Aj j can be expressed in terms of the generators 
j ,j± only and then G(T) = Q(j Q ), F(T) = F(j ). The algebraic morphism j (T) -> W 
and j(T) — > 1 + T s provides the mapping Aj j — > -4. s , a . Since the first two commutators 
of Aj 0j j, eq.(33), do not depend on the function j(T), and the third commutator in eq.(33) 
does not depend on the function jo(T), the closure condition for Aj j is independent of the 
functions jo(T) and j(T) and hence they can be choosen in a convenient way to provide any 
dilation equation in the form of the eigenproblem for j_. The algebraic closure conditions 
for Aj j, together with the definitions of j, jo, require that 

UT 2 )-j (T) = Q(UT 2 )), (36) 

j(T 2 )j(T-i) -j{T-^)j{T) = HUT)). (37) 

Both eqs. (36,37) are non-linear, and in general, difficult to solve analytically. The unique 
solution for the dilation equation are j± since the Casimir operator of this algebra depends 
on T only. Choosing g(T) = j(T) yields the dilation equation in the form 2j_$ = $. In 
the limit of .4.^^ this provides again the Haar two-scale equation. The corresponding scaling 
function belongs to the basis of the representation of j_ with eigenvalue 1/2. The remaining 
arbitrary function j Q can now be selected to obtain the wavelet generating operator j_j $ = 
in agreement with eq.(4). One can formally construct a basis of the representation 
for j_ with the functions \n >= (lnT) n $. The procedure is the following: for a given 
algebra Aj j (therefore given functions Q, J-) solve eqs. (36,37) with respect to the functions 
jo,j± and obtain the corresponding dilation equation in the form 2j„$ = $. A certain 
combination of generators produces the corresponding wavelet equation for ^. Conversely, 
given a scaling/wavelet function, and consequently its dilation equation and g(T), one can 
choose j_ = Dg(T) and then solves the equation Dg(—T~ l ) = j__j with respect to j . With 
jo,j± known one can solve eqs. (36,37) with respect to Q,T to construct the algebra. 

This algebraic approach is in some sense universal. What is modified is the specific 
realization of the non-linear algebra in terms of D and T. Constraints are imposed by the 
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two limiting approaches (two-scale equation and Fourier limit). Loosley speaking, the closure 
of the algebra provides the fixed-point dilation equation and the non-linearity of the algebra 
provides the exponential scaling. 

We also note that Vq = {J2k CfcT fc $},that is, the space generated by all the integer 
translations of $, and any other Vj = D^Vq, is not invariant to the action of j— For 
example, the condition j-DQ = J2k CkT k D§ implies the existence of an operator containing 
D, whose commutator with j_ is a function of T only. From the definition of the algebra 
this is impossible, and consequently this proves the above conjecture. The Vj spaces form, 
by recursion, a basis in L 2 (R). 

In order to prove the uniqueness of the two-scale equation in the general case, we 
again use Proposition 1, for W- — > j_ = Dj^^ 1 ) = DJ2kezJkT h . We have to solve the 
commutation equation [Dj(T~ r ), X)] = for x(T) = EkezX k T k , X = Dx(T). This implies 
that the arbitrary functions rc(£) and j(£) must fulfill the functional condition (£ a generic 
variable) 



Eq.(38) does not carry any restriction with respect to the values of the functions between 
and 1. Since Dj{T) represents a scaling operator we have, according to the dilation equation, 
j(l) = 1, j{— 1) = which implies x(—l) =0. If $ is the corresponding scaling function for 
j_ then X$ is also an eigenfunction of j_, j_(X$) = X$. In order for Dx(T) to satisfy the 
average value property (1) in subsection 2.2, we impose the additional condition x(l) = 1. 
The last condition require of for DX(T), the orthogonality condition (2) of subsection 2.2, 
is equivalent with the condition [1-4] 



By introducing eq.(38) in eq.(39) we obtain X = j_ which provides the trivial indentity 
solution, and hence uniqueness. 



(38) 



x{T)\ 2 + \x{-T)\ 2 = l. 



(39) 
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5 Construction of the scaling and wavelet algebra 



The link between the scaling function and wavelet, and the corresponding scaling algebra is 
supported by the solutions of eqs. (36,37). In the following we present a method for solving 
these equations, from the wavelets towards its algebra. The dilation equation eq. (2) has the 
form of a fixed-point equation. Therefore we must try to find its eigenvectors as limits of some 
functional sequences. The limits of these sequences should be compact supported or rapidly 
decreasing functions for two reasons: to obtain scaling functions with good localization, and 
to provide correct behaviour of the action of the operators f(T). Since the scaling function 
will be expressed as a limit of a sequence, one has to look for operators which commute with 
the limit (are absorbed in the limit). We choose a test function A (x) and a sequence of 
operators f n from the universal covering U(Aj j) of the scaling function algebra introduced 
in the preceding section, such that the limit <3>(x) = lim^oo f n A (x) exists in the weak 
topology [22] and it provides the scaling function. 

In the following we show how to construct the scaling algebra Aj j from the dilation 
equation, written in the form j-(T, D, s)$ = $. We have the following proposition in the 
framework of the algebra Aj j, for s = 1: 



Proposition 2 Let f n {jo) be a functional sequence in U(Aj j), s — 1, A (x) a (test) 
function and <&(x) the scaling function. If the limit $(x) = lim^oo f n (jo)j- Aq(x) exists, 
and lim^oo / "°; , ~ g / s 'f ,)) = I, then j-$(x) = $(ar). 

Jn—1\30) 

Proof: 

From the RHS of the first condition in the hypothesis and from eq.(33) we have 



Jn-ltf0j 

It follows: 

j_$ = fim j_f n (j )A (x) = $, 

q.e.d. It follows from Proposition 2 that for a given dilation equation j_$ = $, (both j_ and 
$ given) we can find out a sequence f n (jo), the operator Q(s, jo), and a test function A (x), 
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such that these objects satisfy the hypothesis of Proposition 2. Then it follows that one can 
construct the scaling algebra, since with j Q and Q(s,j ) found, J-{s, j ) results from the last 
commutator in eq.(33). In general, one starts with j as an arbitrary function of T which 
can be deformed into —id when mapping Aj Q j to JF 01 . As j_ is provided by the two-scale 
equation, j_<& = $, other solutions are forbiden. If lim^^ f n (jo) = f<x>(jo) 7^ const., this 
limit should commute with j_ which is forbidden by Proposition 1. This procedure closes 
the construction of the algebra Aj j. 

In the following we give an algorithm for finding f n (jo) and jo(^) an d illustrate it 
with two examples. Since any mother scaling function is defined by the dilation equation 
(Dh(T)$ = $ or j_$ = $) it is natural to search for solutions by using a recursion algorithm. 
For instance, we can find a test function A (x) and a triplet of operators A, B, C such that 
the limit = B lim^^ A n A (x) exists, and in addition, CB = BA k for a finite positive 
integer k. Then, we have the property C$(x) = §(x). Indeed, C<3> = CB lim^oo A n A = 
B lim^oo A n+k Ao = B lim^oo A n Ao = From the dilation equation we know that C 
should have the form C = Dc(T) with c(T) a function of T. One of the simplest choices is 
to use A = D and B = b(T). Then we have CB = Dc(T)b{T) = c(T 1 / 2 )6(T 1 / 2 ) J D = b(T)A; 
that is, k — 1 and a restriction for the arbitrary functions c(T) and b(T) arises: 



This equation is useful in both directions (algebraic <-> dilation equation), since one can start 
with a given scaling function (given c(T)) and find the operator b(T) involved in the algebra 
and conversely. 

Now consider the algorithm for the Haar scaling function, c(T) = 1 + T _1 . We look 
for solutions of eq.(40) as Laurent series for b(T) = J2kezbkT k . In this case eq.(40) has 
an unique solution, b(T) = const, x (1 — T _1 ). This gives again an unique solution for 
$ = (1 — T _1 ) lim n ^ oo D n A . We stress that Db(T) is exactly the operator which gives the 
Haar wavelet. Further, we can express §Haar in the form = H(x) — H(x — 1), where 

H(x) is the Heaviside distribution. For any sequence of C°° functions S n (x) — > 5(x), we have 
A n (x) = / 5 n (x)dx — > H{x) and we can write 



c(T 1/2 )6(T 1/2 ) = b(T) 



(40) 



®( x )= lim(A n (x)-A n (x-l)). 



(41) 



17 



For example, we can use the sequences 5 n = ^ 2 ^ 2 " 2n or 5 n = co ^^ nx) and A n = aTCta < in ' £ ) 
or A n = ^tanh(2 n x), respectively. The latter example is a soliton-like shape, having good 
localization. These sequences converge to S(x), respectively. In this way we have selected 
subsequences that step in powers of 2. Eq.(41) can be written as: 

= lim (1 - T _1 )D n A (a:), (42) 

n^oo 

We can express this definition in terms of the algebra .4.1,1. By using eq.(25) and the prop- 
erties of the operators D and T, we can write eq.(42) in the form 



$(x) = Jim (1 - T~ )(2W^) n A (x). (43) 
Indeed, from the commutation relation between D and T and the definition of W- we have 
{2W_) n = (D(l + T- l )) n = D n (l + T- 1 + ... + T- 2n+1 ) (44) 

= (i + T ~ 2 ' n + {T- 2 '") 2 + ... + (T~ 2 ~ n ) 2n ~ 1 )D n . 

Hence, we can write (1 — T _1 )D™ = (1 — T~ 2 n )(2WJ) n . Moreover, we can write, by using 
the inverted form for W-(T), the dilation equation in a pure algebraic form: 

= lim (1 - T- 2 ~ n (W ))(2W_) n A (x). (45) 

From this last equation it follows that 2V4 / _$ = $, i.e., the Haar dilation equation. If the 
function Aq(x) is choosen from a class of suitable functions for a wavelet analysis then its 
scaling function <&(x), eq.(45), gives a rapidly convergent wavelet expansion. In order to 
demonstrate this procedure, we show in Fig. 2 scaling functions obtained in this way for 
different values of the parameter s within A s , a - For s = and 1 this functions yields the 
Haar scaling function. 

We give another example for the $ 2 B-scaling function and the \l/ 2 B- wavelet. We 

define 

$ 2 (ar) =4x for < x < 1/2; -Ax + 4 for 1/2 < x < 1, (46) 
and in the rest, or 

*2(x) = Jim - 2T^ 2 + T-^DVViOr), (47) 
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with 

/ x 2a; , n , ln(l + 4a; 2 ) r 2 ^ . . , , 
Aotri(^) = — arctan(2x) - = / A (x)dx (48) 

7T 27T JO 

Following the above algorithm we obtain the corresponding dilation equation 

3- 



$ 2 (x) = ^(l + T- 1 / 2 ) 2 ^^). (49) 



By using the commutation relation dD = 2Dd and the dilation equation eq.(49) we can 
obtain a self-contained equation for the wavelet: 

% ri (x) = D(l + 2T- 1 / 2 + T 2 fm tn . (50) 

The corresponding algebras for Haar and $ 2 B-scaling functions are different. For the B- 
wavelet we have, in A\ ; i 

G(W ) = W - W 2 - l - (51) 

This new algebra fulfills different commutation relations and consequently has a spectrum 
for Wq which is different from that for the Haar scaling algebra. The algorithm introduced 
above can also be used to generate other scaling functions and their corresponding algebras. 
An interesting track is to obtain polygon-like wavelets or, by following the same procedure 
which guided us to obtain the B-scaling function from the Haar scaling function, to obtain 
smoother scaling functions. 



6 The su q {2) limit of Aj 0u 



The quantum group su q (2), or its extensions [21-26,28-31], are associative algebras over C 
generated by three operators, Jo = (Jo)*, J+, an d J_ = («/+)*, satisfying the commutation 
relations 

[Jo, J+] = J+, [J , J-] = - J-, [J+, J-] = [Jol, (52) 
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where [Jo] is a deformation of Jo, that is, a real parameter-dependent function of Jo, holo- 
morphic in the neighbourhood of zero, and approaching 2J in the limit s — > 0. We have 
to find a deforming functional that transforms the Aj j generators into operators satisfying 
the commutation relations eq.(52). For this the first equation in eqs.(33) can be written as 

(jo - G{s, J ))j+ = j+jo, (53) 
and hence, for every entire function we can write 

p(jo ~ G(s,jo))j+ = j+p(jo)- ( 54 ) 
Let us consider the functional equation 

r(£-0(s,O) = r(O-i (55) 

for a given function r(£) with £ generic. If this equation has a solution r(£) that is an entire 
function, then eq.(55) can be written, for = r(£), in the form [T(j ),j + ] = j+, and 
correspondingly [r(j ), j-] = —j-- These equations give the correspondence between the 
wavelet algebras and the q-deformed algebras like su q (2) or other extensions of them. They 
map the structure of Aj j ± , eq.(33), onto the structure described by the first two equations 
in eq.(52). For the third commutator of each of these algebras we use the function r(£), 
which allows eqs.(33) to be reduced to the third equation in eqs.(52) through the mapping: 

Jo = r(j ), J±=j±, (56) 

If T is invertible, the third equation in eqs.(33) in Aj j can be reduced to the third equation 
of eqs.(52) of an algebra defined by the function [J ] with the identification T o = [ 
where []oT means the composition of the two functions, i.e., ([ ] o T)(^) = [r(^)]. In the 
case of su g (2), i.e. [Jo[= [2Jo] g , the function J-{s, jo, s) becomes: 

F(sk) = (57) 

where q = e s , (f>(£) = q~ r ^ has to satisfy the equation 0(£ — G(s,£)) = 90(0- Consequently, 
eqs. (53-55) provide the reduction of su q {2) into Aj j through T. For details of the above 
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technique of mapping and reduction of non-linear algebras one can see examples in [31]. In 
the case of the algebra eqs. (53-55) give the condition 



p(W -l) = W -G{l,W ), 



(58) 



which has the solution 



Jo=p(W ) 



cosh(2- Wo b), 



(59) 



for any arbitrary constant b. By inverting eq.(59) we obtain 



(60) 



Conversely, starting from su q {2) and its deformation [ ] q one can map it into A\,i where for 
the operator F(s,T) = F(s,T(W )) we have 



F = f(T(W )) = f(p-\W )) = [2p-\W Q )\. 



(61) 



A diagram containing these algebraic morphisms is presented in Fig. 3. 



7 Further extensions, comments and conclusions 



In this paper we present a method for connecting classes of non-linear (q-deformed) algebras 
with scaling functions and wavelets. The connection is reciprocal. The algebra Aj 0i j\ s =i 
allows the choice of different dilation equations, i.e., for different functions j(l,T) one can 
obtain different scaling functions. The function jo(l,T) should be choosen in order to give 
the correct wavelet function through the action of joj-. The connection with quantum 
groups is more transparent if we choose for j = [d] s=1 / 2 = W (2, 1/2), eq.(14), and for $, 
the Haar scaling function, j_ = D(l + T _1 ). The action of the q-derivative [<9]i/2 on $ 
gives exactly the Haar wavelet. In this example the operator j is anti-Hermitian and we 
can identify its spectrum (imaginary eigenvalues) from the commutation relations. We also 
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have the relation \l> = (j — G(l,j )&, so that one can express the eigenvectors \a > of J 
in the wavelet basis of the algebra. If \a >= J2j,nCj, n &j,n, we obtain a recursion relation for 
the coefficients of |o > and for the eigenf unctions of j Y^k^j^^k^k = a Cj tP for any j and 
p integer, where are the Taylor coefficients of the function — £ _A <?(— £)■ We note that by 
using the q-derivative instead of a function of T we obtain exactly the same mother wavelet 
(Haar wavelet). 

We note a natural extension of the above developments coming from a realization of 
sp(2, R) in terms of the full algebra of the symmetry of the real line: X = xd (dilations), 
X_ = d (translations), and X + = x 2 d (expansions). This is the maximal finitely generated 
real Lie algebra with generators in the form x n d. By exponentiation the generators X , X_ 
provide D and T and the last generator has the action E a = e aX+ f(x) = f(x/l — ax). 
For a > the action corresponds to a contraction of the function and for — 1 < a < 
it is a dilation. But for a < — 1 and \a\ > 1 we obtain a splitting of the function into 
two Heaviside distributions: E a f(x) = H(x) + H(—x + 1/a). A possible closed non-linear 
algebra containing D, T and E might be an opportunity for the application of wavelets with 
multi-scale properties [18,19]. 

In conclusion, we have found that the operators of dilation and translation (D,T) 
can be combined in such a way as to generate non-linear algebras which depend on certain 
parameters (s,a). We have investigated these algebras from the point of view of quantum 
groups, discussing their unirreps, Casimir operators, and reductions of these algebras to 
other q-deformed algebras, like su q (2) or to the Fourier series generating algebra. It has been 
shown that such algebras provide an appropriate framework for the foundation of wavelet 
analyses and for the obtaining the corresponding scaling functions. We have worked out two 
examples: the Haar and the B-scaling functions. The algorithm provides a general algebraic 
method for finding a specific scaling function. Direct applications of such an approach can 
be found in the theory of finite-difference equations and q-difference equations [21-25,27]. 

This study represents only a first step towards understanding the relation between 
the theory of non-linear algebras (having exponential spectra) and wavelets with their finite- 
difference equations. We conjecture that such an approach to scale invariant structures can 
lead to interesting mathematical constructions and tools for identifying isolated coherent 
structures associated with certain mother wavelets and limiting non-linear algebraic struc- 
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tures. The transition between such self-organized structures can be carried out through the 
modification of the deformation parameter q, with the intermediate domain representing 
"noise" (non-closed algebraic structures). 
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FIGURE CAPTIONS 



1. The dependence on a of three eigenvalues a n of W in A±,i, for n = 2 (the largest 
period), 4 and 6 (the smallest period). The three spectra have self-similar structure 
for a e [0, 1]. 

2. The continuous deformation of the scaling function $ in A s , a as a function of the 
parameter s. For s = 0, 1 this is the Haar scaling function. 

3. A diagrammatic representation of the algebraic maps (arrows) and morphisms (~) of 
Aj oj , A s>a , su q (2), su(2) and F 0jl . 
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